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ABSTRACT 

We present an analysis of '^5 years of Lick Observatory radial velocity measurements targeting a 
uniform sample of 31 intermediate-mass subgiants (1.5 < M*/Mq < 2.0) with the goal of measuring 
the occurrence rate of Jovian planets around (evolved) A-type stars and comparing the distributions of 
their orbital and physical characteristics to those of planets around Sun-like stars. We provide updated 
orbital solutions incorporating new radial velocity measurements for five known planet-hosting stars 
in our sample; uncertainties in the fitted parameters are assessed using a Markov Chain Monte Carlo 
method. The frequency of Jovian planets interior to 3 AU is 26^g%, which is significantly higher 
than the 5-10% frequency observed around solar-mass stars. The median detection threshold for 
our sample includes minimum masses down to {0.2, 0.3, 0.5, 0.6, 1.3} Mjup within {0.1, 0.3, 0.6, 
1.0, 3.0} AU. To compare the properties of planets around intermediate-mass stars to those around 
solar-mass stars we synthesize a population of planets based on the parametric relationship dN oc 
M^P^dlnMdlnP, the observed planet frequency, and the detection limits we derived. We find that 
the values of a and (3 for planets around solar-type stars from Cumming et al. fail to reproduce the 
observed properties of planets in our sample at the 4 a level, even when accounting for the different 
planet occurrence rates. Thus, the properties of planets around A stars are markedly different than 
those around Sun-like stars, suggesting that only a small (^ 50%) increase in stellar mass has a large 
influence on the formation and orbital evolution of planets. 

Subject headings: planetary systems: formation — stars: individual (HD 167042, HD 192699, HD 
210702, K CrB, and 6 Lyn) — techniques: radial velocities 
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1. INTRODUCTION 

Our understanding of planet formation has rapidly im- 
proved over the past 15 years. Prior to the discovery of 
the firs t extrasolar planet orbi ting a solar-type star (51 
Peg b; iMavor fc Ouelod Il995f ). it was widely assumed 
that extrasolar giant planet semimajor axes would mimic 
those of the gas giants in our own solar system, which 
orbit at distances > 5 AU. In the years that followed 
it became apparent that an in situ formation model was 
not universally applicable because radial velocity surveys 
were finding Jovian planets in abundance well inside the 
canonical ice line. 

Over 350 planets have now been discovered, 282 of 
which reside around stars within 200 pc.^ Sufficiently 
large samples are available for the statistical proper- 

^ Based on observations obtained at the Lick Observatory, which 
is operated by the University of California 

^ Institute for Astronomy, University of Hawai'i, 2680 Woodlawn 
Drive, Honolulu, HI 96822, USA; bpbowlerOif a.hawail . edu 

^ Current address: Department of Astrophysics, California In- 
stitute of Technology, MC 249-17, Pasadena, CA 91125, USA 

* Department of Astronomy, MS 3411, University of California, 
Berkeley, CA 94720-3411, USA 

^ Center of Excellence in Information Systems, Tennessee State 
University, 3500 John A. Merritt Blvd., Box 9501, Nashville, TN 
37209, USA 

® Department of Physics and Astronomy, San Francisco State 
University, San Francisco, CA 94132, USA 

''' ZAH-Landessternwarte, Konigstuhl 12, 69117 Heidelberg, Ger- 
many 

8 UCO/Lick Observatory, Santa Cru z, CA 95064 
® As of August 2009; see [http: / /exoplanets.org| 



ties of exoplanets to reveal themselves, providing infor- 
mation about the planet formation and migration pro- 
cesses. For solar-type stars (F, G, and K dwarfs), Jo- 
vian planets fall into two rough populations: "hot plan- 
ets" with a < 0.1 AU (Fischer 200§) and those that 
orbit beyond ^ 1 AU (Figure [T]) . These observations 
were explained a posteri ori in terms of orbital migration 

Paoaloizou et al.l 120071) and planet-planet scatterin; 

Nagasawa et al.ll2008l : iMarchi et al.ll2009l : iFord fc Rash 



like stars (e.g.. \tiu 
iMarcv et al.l l2008t 



2008') , with the dearth of planets with periods between 
~ 10-100 days (the "period valley") possibly result- 
ing from differen t ial mass-dependent orbital migra tion 
(|Udrv et al.ll200l IBurkert fc Idal[2007l : ICurriell2009l) . 
While much is known about planets around Sun- 
.IButler et all 120061: lUdrv fc SantosI 120071 : 
' JohnsonI 120091 : I Wright et al.ll2009fl . 
comparatively little is known about planets around 
intermediate-mass (IM) stars with > 1.5 Mq. Main 
sequence A- and F-type stars are problematic Doppler 
targets because of their high jitt er levels and rotation ally- 
broadened absorption features (jGalland et al.ll2005[ ). On 
the other hand, evolved IM stars have lower jitter 
levels as well as narrower and more numerous ab- 
sorption lines resulting from their slow rotation and 
cool photospheres. As a consequence, nearly all ra- 
dial velocity surveys of IM stars are targeting evolved 
G- or K-type subgiants and gian ts (Frink et al. 20Q3, 
Setiawan et al. 2003a, Sato et al. '2005b: Johnson et aT] 
2006, Nicdziclski ct_al. 2007, DoUingcr et al. 200'^ 



Lovis fc Mavoiil200'7l . iLiu et all 120091) . although at least 
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Fig. 1. — The semimajor axis distribution of exoplanets within 
200 pc. The data are compiled from the literature and are main- 
tained by the California & Carnegie Planet Search team. Planet- 
hosting stars with masses less than 1.5 Mq are shown as open 
squares, while planet-hosting stars with higher masses are shown as 
filled diamonds. For solar-type stars, planets can be grouped into 
two populations: "hot planets" located at distances ^0.1 AU and 
longer period planets located at distance > 1 AU. Between these 
two groups lies a relative dearth of planets in a region known as the 
"period valley." No planets have been discovered with semimajor 
axes < 0.6 AU for stellar masses > 1.5 Mq, creating a "planet 
desert" in that region. Planets with semimajor axes > 5 AU are 
limited by the current baseline of radial velocity observations. 

one survey is targeting their main sequenc e progenitors 
(|Galland e t al. 2005; Lagra nge et al.ll2009aD . 

Recently .Johnson et al.. (|2007b[ ) showed that planets 
orbiting evolved A-type stars have large semimajor axes 
compared to planets around solar-type stars. This trend 
has become even stronger with the discover y of more 
planets orbiting IM stars (jSato et al.l I2008af) . Specifi- 
cally, over 20 planets have been discovered around stars 
with minimum masses > 1.5 Mq but none have semi- 
major axes < 0.6 AU. This "planet desert" is shown in 
Figure [T] and a summary of known planets around IM 
stars is presented Table [TJ 

There are several explanations for this observational 
result. One possibility is that the swollen radii of evolved 
stars have engulfed or tidally disrupted short-period 
planets. The current census of IM exoplanet-host stars 
mostly consists of red clump giants (core helium and 
hydrogen shell burning post-RGB stars) and subgiants 
(stars with contracting inert helium cores and hydrogen- 
burning shells). Planet-hosting stars in the clump giant 
phase have radii between - 8-14 Rq (- 0.03-0.06 AU) 
so any planets orbiting inside ~ 0.1 AU will likely have 
been engulfed by an expanding radius during post-main 
sequence stellar evolution. However, a tidal torque from 
an expanding stellar su rface can a l so deca y the orbits 
of short period planets. iSato et all (|2008aD numerically 
trace the semimajor axis evolution of short-period plan- 
ets around evolving RGB stars and show that, in their 
past, IM clump giants may have engulfed or disrupted the 
orbits of planets out to ~ 0.5 AU. With radii between ~ 
2-5 Rq {r^ 0.01-0.02 AU), planets orbiting IM subgiants 
are the least affected by stellar evolution. Stellar evolu- 
tion may therefore explain the lack of short-period plan- 
ets around IM clump giants, but the same result for IM 
subgiants suggests that the observed trend is not due to 



post-main sequence engulfment. 

The observed semimajor axis distribution can also be 
explained as a result of inward orbital migration com- 
bined with mass-dependent disk dispersal lifetime. The 
rocky progenitors of Jovian planets (~10 Af® cores) can 
form at distances > 8 AU for s tellar m asses b e tween 
~ 1.5-3 Mq (IKennedv fc KenvonI l200l . ICurrid (l2009l ) 
performed Monte Carlo simulations of Jovian planet for- 
mation and migration around IM stars and, using sim- 
ple stellar mass-dependent gas disk lifetime relations 
(iKe nnedv fc Kenyon 2009) and Type II migration mod- 
els (jida fc Linil2f)0l . was able to successfully reproduce 
the observed dearth of short-period planets. In this sce- 
nario inward migration is halted once rapid disk dispersal 
occurs, stranding migrating planets at semimajor axes 
that depend on stellar mass. 

An alternative explanation was offered bv lKretke et al.l 
(|2008f l. They found that the protoplanetary disks of 
young IM stars will develop a maximum surface density 
at ~ 1 AU as a result of magnetrotational instability of 
the inner disk, leading to a trapping and accumulation of 
solids that can then grow to form rocky cores and Jovian 
planets. This formation scenario provides a mechanism 
for in situ formation of Jovian planets interior to the ice 
line, which is located near 3 A U at 10 Myr for a 2 AIq 
star (jKennedv fc Kenvon|[200l . 

Models of Jovian planet formation around IM stars 
make few quantitative predictions that can be obser- 
vationally tested. Simple disk depletion plus migration 
models for IM planet hosts (1.5-3.0 Mq) bv lCurriel (pOOOl ) 
predict occurrence rates for Jovian planets with semi- 
major axes < 0.5 AU to be < 1.5% and for Jovian 
planets with semini a ior ax es > 0.5 AU to be > 7.5%. 
IKennedv fc KenyonI ^00§) use a semianalytic model of 
protoplanetary disk evolution to study snow line loca- 
tions and planet formation rates around stars of vary- 
ing masses. For stellar masses between 1.5-2.0 Mq, 
their models predict that Jovia n planet occu r rence rates 
reach frequencies of '^10-15%. iKretke et al.l (|2008D sug- 
gest that multiple planetary systems may form more ef- 
ficeintly around IM stars compared to other stellar mass 
regimes. 

The semimajor axis distribution of planets with mini- 
mum masses between ~ 2-10 Mjup orbiting IM stars is 
beginning to be better constrained by observations, but 
little is known about planets with masses < 1.5 Mjup. 
This poor understanding is a direct result of the dearth 
of low-mass planets currently known, with only one hav- 
ing a minimum mass below 1.5 Mjup (HD 175541b with 
Mpsim = 0.70 Mjup). It is unclear, however, whether 
this scarcity is a result of a detection bias caused by 
higher jitter levels in IM stars or whether it reflects an 
intrinsic shortage of low-mass planets. Unfortunately, 
models of planet formation in this stellar mass regime 
have made few predictions about low-mass planetary 
companions. For solar-type stars, iMavor et al.l (|2009[ ) 
estimate the frequency of Neptune-mass planets with pe- 
riods < 50 days to be at least 30%. This raises the excit- 
ing possibility that low-mass planets could be abundant 
around IM stars, especially in light of recent studies sug- 
gesting that the fr equency of Jovian-ma ss planets scales 
with stellar mass pohnson et al.|[2007af ). Testing these 
theories requires an understanding of the detection limits 
of radial velocity surveys. 
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TABLE 1 

Planetary-Mass Companions to Evolved Intermediate-Mass Stars with Af* > 1.5 Mq 
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Star 




SpT 




V, ' J up / 


(AU) 


e 


Ref 


HD 13189 


2-6 


K2 H 




8-20 


1.5-2.2 


0.27 (0.06) 


1, 2 


e Tau 


2.7 (0.1) 


KO HI 


13.7 (0.6) 


7.6 (0.2) 


1.93 (0.03) 


0.151 (0.023) 


3 


NGC 2423 No. 3 


2.4 (0.2) 






10.6 


2.10 


0.21 (0.07) 


4 


81 Get 


2.4 (2.0-2.5) 


G5 HI: 


11 (10-13) 


5.3 


2.5 


0.206 (0.029) 


5 


HD 104985 


2.3 


G9 HI 


11 


8.3 


0.95 


0.090 (0.009) 


6, 7 


18 Del 


2.3 


G6 III 


8.5 


10.3 


2.6 


0.08 (0.01) 


6 


HD 17092 


2.3 (0.3) 


KO III 


10.9 (2.8) 


4.6 (0.3) 


1.29 (0.05) 


0.166 (0.052) 


8 


€ Aql 


2.2 


KO III 


12 


2.8 


0.68 


0.0 (fixed) 


6 


14 And 


2.2 (2.0-2.3) 


KO HI 


11 (10-12) 


4.8 


0.83 


0.0 (fixed) 


5 


HD 81688 


2.1 


KOIII-IV 


13 


2.7 


0.81 


0.0 (fixed) 


6 


HD 173416 


2.0 (0.3) 


G8 HI 


13.5 (0.9) 


2.7 (0.3) 


1.16 (0.06) 


0.21 (0.04) 


9 


HD 11977 


1.91 (0.21) 


G5 HI 


10.09 (0.32) 


6.54 


1.93 


0.40 (0.07) 


10, 11 


HD 102272 


1.9 (0.3) 


KO HI 


10.1 (4.6) 


5.9 (0.2) 


0.614 (0.001) 


0.05 (0.04) 


12 








2.6 (0.4) 


1.57 (0.05) 


0.68 (0.06) 


12 


13 Gem 


1.86, 1.7 (0.4) 


KO III 


8.8 (0.1) 


2.9 (0.3) 


1.69 (0.03) 


0.06 (0.04) 


13, 14, 15, 16 


HD 89744 


1.86 (0.18) 


F7 IV 


2.08 (0.06) 


7.2 


0.88 


0.70 (0.02) 


17, 18 


HD 210702 


1.85 (0.13) 


Kl IV 


4.45 (0.07) 


1.97 (0.11.0.18) 


1.20 (0.02,0.03) 


0.036 (<0.106) 


19, 20 


K GrB 


1.84 (0.13) 


KO IV 


4.71 (0.08) 


2.01 (0.11,0.17) 


2.80 (0.07,0.08) 


0.044 (<0.123) 


21, 20 


6 Lyn 


1.82 (0.13) 


KO IV 


5.2 (4.9-5.6) 


2.21 (0.11,0.16) 


2.18 (0.05,0.06) 


0.059 (<0.125) 


20, 5 


HD 167042 


1.72 (0.12) 


Kl IV 


4.30 (0.07) 


1.70 (0.09,0.12) 


1.32 (0.03,0.04) 


0.089 (0.028,0.065) 


21, 20, 5 


HD 192699 


1.69 (0.12) 


G8 IV 


3.90 (0.06) 


2.40 (0.15,0.21) 


1.15 (0.02,0.03) 


0.129 (0.029,0.060) 


19, 20 


HD 175541 


1.65 (0.12) 


G8 IV 


3.80 (0.09) 


0.70 (0.06,0.08) 


1.03 (0.02,0.03) 


0.083 (<0.283) 


19 


HD 5319 


1.59 (0.18) 


G5 IV 


3.26 (0.50,0.41) 


1.94 


1.75 


0.12 (0.08) 


22 


References. — (1 


iHatzcs ct aL ( 


200,'5); (2) Schulcr ct aL (2005); 


(3) Sato ct aL (2007): (4) Lovis & Mayor (2007); (5) Sato ct aL (2008b); (6) 



ISato et al] (|2008aD :_f 7) fSato ct aL (2003); (8) Nicdziclski ct aL 12007); (9) Liu ct aL (200 9); (10) Sctiaw an ct aLJ2005); (11) da Silva ct aL ( 200(^') ; 
(12) Nicdziclski ct^al.' ( 20091 ) ; (13) J latzcs ct al. (200d); ('14 ) INordgren et al.l (120011) ; (15) |Prieto fc LambertI (Il999j); fl 6 )iReffert c t al. (200 Q); (17) 
iKorzcimik ct al. (2000); ( 18) IValent i & Fischer (20(^; (19) I Johnson et al.l l l2007bl) ; (20) this work: (21) ljQhnson et al.l (120081) ; (22) Robins on et al'l 
1200711 



No te. — HD 47536 b is omitted because the most likely stellar mass is between 1.0-1.5 (ISetiawan et al.ll2003al) . 7 Cep A b (IHatzes et al.l 

1200311 is omitted because the stellar mass estimate was revised to 1.4 Mq lINeuhauser et al.l2007^ . The minimum mass of the companion to 11 Comae 
is 19.4 Afjup JLiu c t al. 2008) and is therefore probably a brown dwarf. References are for planet discoveries and host-star physical properties. 
Additionally, HD 90043 and HD 200964 arc intcrmediatc-mass planet-hosting stars, but do not yet have unique orbit solutions and so were omitted 
from this table. 



The goal of this study is twofold: to derive the oc- 
currence rate of Jovian planets around IM stars and to 
characterize the distributions of planet periods (P) and 
minimum masses (Afpsim). In addition we take the op- 
portunity to update the orbit solutions for known planet- 
hosting stars in our sample using new radial velocity mea- 
surements. To address the aforementioned questions we 
make use of a uniform sample of 31 IM subgiants taken 
from an ongoing radial velocity survey at Lick Observa- 
tory (Johnson et al. 2006; Peek et al. 2009). Our obser- 
vations span ~ 5 years and sample semimajor axes out to 
several AU, enabling detailed comparisons with previous 
surveys targeting solar-type stars. 

In §2 we describe our radial velocity measurements and 
define our sample selection. Updated Keplerian orbits for 
five previously-known planetary systems are presented in 
§3 and in §4 we describe time-series photometric observa- 
tions of these systems. In §5 we derive detection limits for 
stars in our sample. We discuss the frequency of Jovian- 
mass planets in §6, and in §7 we compare the mass-period 
distribution for solar-type stars to the results from our 
sample of IM subgiants. Finally, we discuss the implica- 
tions of our work in §8. 

2. RADIAL VELOGITY OBSERVATIONS 

Our sample of subgiants is derived from a larger planet- 
hunting program at Lick Observatory tar geting evolved 
stars (the Lick Subgian t Planet Search: iJohnson et al.l 
I2006t iPeek et al.ll2009[ ). The original survey consists 
of 159 stars selected on the basis of Hipparcos B — V 



colors and My absolute magnitudes such that they lie 
one magnitude above the main sequence, have masses > 
1.2 Mq, and avoid red giant branch (RGB), clump gi- 
ant, and Cepheid variable regions of the HR diagram. 
Known spectroscopic and proper motion binaries within 
^ 2" are excluded. To obtain a more uniform popula- 
tion for this study we selected a sample of stars from the 
original survey based on the following criteria: 2 < My 
< 3, 0.8 < B - V < 1.0, Af, > 1.5 Mq, and at least 
four radial velocity measurements (the average number 
of measurements in our final sample is 26).^° The magni- 
tude cuts exclude clump giants and solar-mass subgiants, 
while the color cuts exclude RGB stars and yellow strag- 
glers, the latter of which have high jitter levels and large 
binary frequencies. Altogether our sample includes 31 
intermediate-mass subgiants (Figure [21). 

Radial velocity measurements were obtained at Lick 
Observatory's 3 m Shane telescope and 0.6 m Coude 
Auxiliary Telescope with the Hamilton spectrometer 
(Vod (19871 : R - 50,000 at A = 5500 A). Doppler shifts 
are measured from each spectrum u sing the iodin e cell 
method described in detail by Butler et al.l (|1996D and 
summarized as follows. A temperature-controlled Pyrex 
cell containing gaseous iodine is placed at the entrance 
slit of the spectrometer. The dense set of narrow molec- 
ular lines imprinted on each stellar spectrum from 5000 
to 6000 A provides a robust wavelength scale for each ob- 

Only one star, HD 33066, has < 4 radial velocity measure- 
ments, but it is a previously-known spectroscopic binary. 
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Fig. 2. — Color-magnitude diagram showing the location of 
our subsample of subgiants monitored at Lick Observatory. The 
thick line shows the locus of Hipparcos main sequence stars and 
the thin lines displ ay the solar metallicity evolutionary tracks of 
IGirardi et al.l l|2002l ) . Our uniform sample of 31 subgiants is a sub- 
set of the larger Lick Subgiant Planet Search program. 
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Fig. 3.— Updated Keplerian orbit for HD 167042. The gray 
dashed line shows the best-fitting orbital solution. The lower panel 
displays the residuals after subtracting off the model. There is 
evidence for an additional outer companion based on a linear trend 
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of 2.14^0 ?9 



servation, as well as information about the shape of the 
spect rometer's instrumental response (|Marcv fc Butleii 
119921 ). Additional details abou t the observations can be 
found in I Johnson et al.ll2007bl 

3. KEPLERIAN ORBITS AND UNCERTAINTIES 

We derive Keplerian orbital solutions for the five 
previously-known planets in our sample using RVLIN, 
an efficient orbit-fitting routine written in IDL and 
based on a partial linea rization of Kepler's equations 
(| Wright fc Howard! 120091 ). The radial velocity measure- 
ments were first combined into two-hour bins and the 
total measurement errors were estimated by adding 5 m 
s""'^ to the internal measurement errors, in quadrature, to 
account for the moderate jitter level of typical subgiant s 
(|Fischer et al.l[200l IWrigliH[2005t IJohnson et alll2007bl) . 
We combine our observations of HD 167042 an d 6 Lyn 
(HD 45410) with those made available by iSato et all 



Fig. 4.— Updated Keplerian orbit for HD 192699. The gray 
dashed line shows the best-fitting orbital solution. The lower panel 
displays the residuals after subtracting off the model. 
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Fig. 5.— Updated Keplerian orbit for HD 210702. The gray 
dashed line shows the best-fitting orbital solution. The lower panel 
displays the residuals after subtracting off the model. 
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Fig. 6.— Updated Keplerian orbit for k CrB (HD 142091). The 
gray dashed line shows the best-fitting orbital solution. The lower 
displays the residuals after subtracting off the model. 
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Fig. 7.— Updated Keplerian orbit for 6 Lyn (HD 45410). The 
gray dashed Une shows the best-fitting orbital solution. The lower 
panel displays the residuals after subtracting off the model. 



((2008b) to create a larger data set for improved orbital 
phase coverage. The data sets were merged by simulta- 
neously fitting for the constant offset between the two 
observatories. 

Long-period companions will manifest as constant ac- 
celerations in the radial velocity measurements. To first 
order this drift can be modeled as a linear velocity trend 
for very long-period companions for which the period is 
much greater than the observing time baseline. When 
the slope has an amplitude comparable to the noise it 
is often difficult to determine whether a trend should 
be included in an orbit solution. The goodness of fit 
statistic (typically the value) might be reduced with 
the orbit-plus-trend solution, but, by that criterion for 
model selection, a more complicated model with multi- 
ple planets could be fit to a radial velocity set for which 
there is clearly insufficient evidence for additional com- 
panions. More generally, how does one decide whether 
including more parameters in a model is justified? For- 
tunately this problem has received a great deal of atten- 
tion in statistics and is the subject of Bayesian Model 
Comparison. A commonly used model s election tool is 
the Bayesian Information Criterion (BIC: ISchwardfTOTSl : 
iriddle'2004; Liddlc 2007), which rewards better-fitting 
models but penalizes models that are overly complex: 



BIC = -21n£n 



fclniV, 



(1) 



where C^ax is the maximum likelihood for a particular 
model with k free parameters and N data points. In 
general the model with the smaller BIC value is preferred. 

We use BIC values for orbit solutions with and with- 
out a velocity trend to decide whether to include a lin- 
ear slope in the model. The only planet-hosting star 
for which a trend produced a lower BIC value was HD 
167042. We take this as evidence for a long-period com- 
panion in that system. The best-fitting orbits for the five 
known planet-hosting subgiants are shown in Figures [3]- 
[7l and the updated orbital parameters are listed in Table 

A popular frequentist approach to model selection is to derive 
false alarm proba bilities for a more comp licated model using of the 



r bilities tor a more comp licated 
F-test (see, e.g., ICumming et al.l 119991) . We 
Bayesian approach in this work. 



choose to follow a 



1 . 

Minimum masses and semimajor axes are calculated 
using the analytic approximation for Mpsini as a func- 
tion of the observed orbital paramete rs and Kepler's 
third law (i.e.. Equations 2 and 3 in Cumming et al.l 
1 19991 ). We use newly de termined stellar masses following 
the method described in I Johnson et al.l ((2b07 bl based on 
new [Fe/H] values (see Table [Ij, which are derived from 
updated iodine-free template spectra follo\v i ng th e iter- 
ative scheme from Figure 1 of lValenti etall (|2009h . 

We use a Markov-Chain Monte Carlo (MCMC) algo- 
rithm following the description in iFordI (2005) to esti- 
mate our orbital parameter uncertainties (see also iFordI 
[20061 iFord fc GregOTvl [200l . MCMC is a Bayesian in- 
ference technique that uses the data to explore the shape 
of the likelihood function for each parameter of an input 
model. This method of error analysis has been shown 
to be more efficient and accu rate than b ootstrap Monte 
Carlo resampling techniques (jFordl 120051 ). Every step in 
the chain represents a variation of the period (P) , eccen- 
tricity, (e), time of periastron passage (Tp), argument 
of periastron (uj), velocity semi-amplitude (K), and con- 
stant velocity offset. A velocity trend parameter (dw/dt) 
is allowed to vary for HD 167042. Additionally, for HD 
167042 and 6 Lyn, a constant radial velocity offset be- 
tween Lick Observ atory measurements and those from 
iSato et aTl (|2008bf) is treated as a free parameter to ac- 
count for the systematic offset between observations from 
different telescopes. We allow one random parameter to 
be altered at each step by drawing a new value from 
a Gaussian transition distribution. If the resulting 

value for the trial orbit is less than the previous 
value, then the trial orbital parameters are retained. If 
not, then the probability of adopting the new value is 
equal to the ratio of the probabilities from the previ- 
ous and trial ste ps (the Metrop olis-Hastings algorithm; 
iMetroDofis et al.lll953t lHastingsiri970 '). If the trial is re- 
jected then the parameters from the previous step are 
adopted. 

We altered the standard deviations of the Gaussian 
transition functions so that the acceptance rates were be- 
tween 20% and 40% to maximize convergence efficiency. 
For each planet-star system, the initial parameters were 
chosen from the best-fitting orbit solutions and each 
chain was run for 1-3 x 10^ steps depending on the rate of 
convergence. The initial 10% of each chain was excluded 
from the final estimation of parameter uncertainties. The 
parameters were tested for convergence by running five 
shorter chains with 10^ steps. We verified that conver- 
gence wa s reached by ensuring t hat the Gelman-Rubin 
statistic C Gelman fc R{Ibinl[l99a : TCowles fc Carlinlll996fl 
for each parameter was near unity (typically < 1.02) and 
that a visual inspection of the history plots suggested 
stability. 

The results of the MCMC analysis are posterior prob- 
ability density functions (pdfs) for each orbital param- 
eter used in the model. An example of the results for 
HD 167042 is presented in Figure [51 For each parameter 
we compute the range encompassing 68.3% of the data 
about the median of the distribution. These values rep- 
resent approximate upper and lower one-sigma limits for 



Similar to, but not to be confused with, 
IIGeman fc Geman|[T98l : IPress et al]|2007l) . 



Gibbs sampling 
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TABLE 2 
Updated Orbital Parameters 



cir a,in G t e r 






upl 9107091-1 




C T U 

6 Lyn b 


P (days) 


420.77 (3.48, 3.11) 


345.53 (1.77, 1.63) 


354.29 (2.31, 2.13) 


1261.94 (28.91, 23.97) 


874.774 (16.27, 8.47) 


Tp (JD - 2,450,000) 


4230.1 (40.6, 41.6) 


4036.6 (21.0, 22.1) 


4142.6 (78.3, 100.1) 


3909.2 (332.5, 260.9) 


4024.5 (180.2, 130.9) 




0.089 (0.028, 0.065) 


0.129 (0.029, 0.060) 


0.036 (<0.106) 


0.044 (<0.123) 


0.059 (<0.125) 


I n 


85.7 (35.4, 35.8) 


41.3 (23.3, 25.9) 


282.3 (82.3, 98.0) 


148.4 (83.2, 94.2) 


314.9 (75.5, 52.6) 


K (m s-l) 


32.16 (1.32, 1.32) 


49.3 (2.6, 3.2) 


37.45 (1.90, 2.49) 


25.17 (1.12, 1.55) 


31.53 (1.12, 1.32) 


Mp Sim (Mjup) 


1.70 (0.09, 0.12) 


2.40 (0.15, 0.21) 


1.97 (0.11, 0.18) 


2.01 (0.11, 0.17) 


2.21 (0.11, 0.16) 


a (AU) 


1.32 (0.03, 0.04) 


1.15 (0.02, 0.03) 


1.20 (0.02, 0.03) 


2.80 (0.07, 0.08) 


2.18 (0.05, 0.06) 




1.17 


1.82 


1.49 


1.01 


1.43 


BIG (trend) 


99.84 


104.97 


90.19 


82.51 


110.25 


BIG (no trend) 


101.88 


101.72 


89.18 


80.15 


106.31 


dv/dt (m s""*" yr^-') 


2.14 (0.99, 0.73) 


0.0 (fixed) 


0.0 (fixed) 


0.0 (fixed) 


0.0 (fixed) 



Note. — Updated orbital parameters for planets orbiting intermediate-mass stars in our sample. In parentheses the upper and lower boundaries 
of the regions encompassing 68.3% of the posterior probability density distributions about the median are listed. For low eccentricities, the peak 
of the probability density distribution may not be encompassed by the 68.3% region about the median. In those cases an upper limit encompassing 
95.4% of the data is quoted. 
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P (Days) 
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Tp (JD - 2,450,000) 
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ci) (deg) K (m s"') 



2 4 

di'/dt (m s"' yr') 



Fig. 8. — Markov Ghain Monte Garlo posterior probability den- 
sity distributions for the orbital parameters of HD 167042b. Me- 
dian values of the distributions are indicated with dotted lines while 
the results from the orbit-fitting routine RVLIN are indicated with 
solid lines. Dark and medium gray regions display the area en- 
compassing 68.3% and 95.4% of the data about the median (cor- 
responding to 1 (T and 2 a areas of a Gaussian distribution). 

Gaussian-like distributions and are included in parenthe- 
ses next to each parameter in Table [2] For eccentricities 
near zero, the 68.3% range of data about the median is 
a poor estimator of the most likely values as it fails to 
encompass the peak of the distribution. When that oc- 
curs we quote upper limit eccentricity values below which 
encompass 95.4% of the data. We note that the results 
from using the BIG value as a way to discriminate be- 
tween adopting models with or without velocity trends 
yields the same results as using the posterior pdf of the 
velocity trend from MGMG as an indicator. 

We use Lomb-S cargle (LS) periodograms (|Lombl 
119761 : iScargld [1982D to search for additional short pe- 
riod companions after subtracting the best-fitting or- 
bit from the observations. A false alarm probability 
(FAP) is compute d follo wing the method described in 
iHorne fc Baliunad (|1986f l. which gives the probability 
that a peak will occur by chance assuming the data are 
pure noise. The results are shown in Figure (S] No peri- 
odicities in the residuals are identified with FAP of 0.1% 
or less. 



HD 167042 Residuals 



FAP = 0.1% 
" FApVl°/o~ ' 



HD 392699 Residuals 



15: HD 210701 Residuals 



K CrB Residuals 



6 Lyn Residuals 



100 1000 



10 100 1000 

P (Days) 

Fig. 9. — Lomb-Scargle periodograms for the residuals of five 
planet-hosting subgiants. False alarm probabilities (FAP) of 1% 
and 0.1% are shown as dotted and dashed lines, respectively. No 
significant periodic signals were detected. 



3.1. The HD 8375 SBl System 

Several stars in our sample exhibit very high radial 
velocity scatter indicative of stellar companions orbit- 
ing at small semimajor axes. HD 8375 is one such 
system and, with 28 observations, is well-enough sam- 
pled to derive an accurate orbit solution and poste- 
rior pdfs using RVLIN and MGMG ( Figure ITUl). HP 
8375 is a known spectroscopic binary fBeavers & Fitted 
1986t [delledciros & Mayoil Il999t ISnowden & Youn^ 
20051 ) but to our knowledge no orbit has been published. 
The amplitude and period of the companion suggest that 
we cannot use the negligible companion mass approxi- 
mation to derive M2sinz. Instead we compute the mass 
function fi{m) and minimum separation aisim, which 
arc listed along with the orbital parameters of HD 8375B 
in Table O We find that a companion with a long-term 
velocity trend best-fits the data based on BIG values with 
and without a trend. 

Following I Winn et al.l (|2009bf ). we can estimate the 
minimum mass required to accelerate the HD 8375 sys- 
tem at the observed rate as a function of semimajor axis: 
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Fig. 10. — Lick Observatory radial velocities and best-fitting Ke- 
plerian orbital solution for HD 8375. The best-fitting orbit includes 
the inner companion HD 8375B with a linear velocity trend {dv/dt 
= 52.9^j g m s~^ yr~^), indicative of an additional long-period 
companion. 



TABLE 3 
Orbital Parameters for HD 8375B 



Parameter 



HD 8375B 



P (days) 

Tp (JD - 2,450,000) 

e 

K (m s-i) 
/i(m) (10-3 
aisini (10 ^ 



Mq) 

AU) 



BIG (trend) 
BIG (no trend) 
dv/dt (m s" yr" 



83.9408 (0.0016, 0.0015) 

4039.3 (0.5,0.4) 
0.0179 (0.0004, 0.0004) 
321.83 (1.74,1.97) 

4939.2 (2.6,2.5) 
1.050 (0.002, 0.002) 
3.813 (0.002, 0.002) 
5.69 
119.3 
928.8 
52.9 (1.8,1.9) 



Note. — See Tablc[2]for details about the orbital 
parameters and uncertainties. 



, r ■ . dv X 
Aipsmi ^ — — sec(27rTj, 



(2) 



where Mp is the minimum planet mass, i is the inch- 
nation, dv/dt is the acceleration, G is the gravitational 
constant, and t is the ratio of the time baseline of the ob- 
servations divided by the period of the companion. For 
large periods relative to t he time base line, sec(27rr) ^ 1, 
yielding Equation 1 from I Winn et al. (2009b). For small 
values of r, the measured value of dv/dt gives 



In addition to the radial velocities from Lick Observa- 
tory, we used the T3 0.4 m and the T12 0.8 m Auto- 
mated Photometric Telescopes (APTs) at Fairborn Ob- 
servatory to gather time-series photometry of the five sys- 
tems with orbital updates in Table [21 The observations 
of these five systems were conducted at various epochs 
between 1993 April and 2009 June. The APTs can de- 
tect short-term, low-amplitude brightness variations in 
cool stars due to rotational modulation of the visibil- 
ity of s urface features, such as starspots and plages (e.g., 
iHenry e t al. 1995,), and can also detect longer-ter m varia- 
tions associated with s tellar magnetic cycles fe.g.. lHenrvl 
119991 iHall et all [20091 ). Photometry of planetary candi- 
date host stars helps to establish whether observed low- 
amplitude radial velocity variations are caused b y stel- 
lar activ ity or planetar y-reflex mot ion (e.g . , Hcnr v et alJ 
I2000al) . iQueloz et all (pOOll ) and iPaulson et al.l " (j200l ) 
have found periodic radial velocity variations in solar- 
type stars that are caused by starspots. The APT ob- 
servations are also useful to search for possible transits 
of the planetary companions (e.g., iHenrv et al.l l2000bl : 
iSato et al.l[2005al : iWinn erani2009a( K 

The T12 APT is equipped with a two-channel precision 
photometer employing two EMI 9124QB bi-alkali photo- 
multiplier tubes (PMTs) to make simultaneous measure- 
ments of a star in Stromgren h and y passbands. The 
T4 APT has a single-channel precision photometer that 
uses an EMI 9124QB PMT to measure a star sequen- 
tially through h and y filters. The APTs observe each 
target star (star D) in a quartet with three ostensibly 
constant comparison stars (stars A, B, and C). We com- 
pute h and y differential magnitudes for each of the six 
combinations of the four stars: D — A, D — B, D — C, 
C — A, C — B, and B — A. Wc then correct the Stromgren 
b and y differential magnitudes for differential extinction 
with nightly extinction coefhcients and transform them 
to the Stromgren system with yearly mean transforma- 
tion coefficients. Finally, we combine the Stromgren b 
and y differential magnitudes into a single {b + y)/2 pass- 
band to improve the precision of the observations. Henr^ 
119991 presents a detailed description of the automated 
telescopes and photometers, observing techniques, and 
data reduction procedures needed for long-term, high- 
precision photometry. 

We use the two most constant comparison stars in each 
quartet to compute the quantity aj, the target star's 
absolute night-to-night variability, statistically corrected 
for any unrecognized intrinsic variability in the two com- 
parison stars and also corrected for the measurement un- 
certainty of the diff erential magnitu des. We follow the 
method described in IHall et al.ll2009l and compute cr/ as 



Mpsini ^ 0.295(a/l AV)^Mj^^ 



(3) 



Values of Mpsini and a below this curve are not permit- 
ted; masses are not high enough at large separations to 
induce the observed acceleration. Lowering the value of 
r shifts the curve towards smaller semimajor axes and 
higher masses. If the third body is a planet then it is 
in a circumbinary orbit, otherwise it is a star or brown 
dwarf in a hierarchical configuration. 

4. PHOTOMETRY OF THE FIVE KNOWN PLANETARY 
SYSTEMS 



(4) 



where cr, is the standard deviation of the target star's 
differential magnitudes computed as the mean of the two 
best D~A, D-'B,D-C time series (e.g., D - {A + 
B)/2), (Tc is the standard deviation of the differential 
magnitudes of the two best comparison stars (e.g., B — 
A), and e is the measurement precision of an individual 
differential magnitude. 

We estimate e by examining the standard deviations 
of the comparison star differential magnitudes (C — A, 
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Detection Limits: Method 
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Fig. 11. — Method for computing detection limits. For each pe- 
riod we perform a Monte Carlo simulation of radial velocity mea- 
surements at the exact dates of the observations beginning with 
a large velocity semiamplitude and gradually lowering it until the 
simulated distribution of rms values is consistent with the observed 
value for that star. For each period and velocity semiamplitude we 
test for consistency at some confidence level; in this study we use 
95.4%. If the observed value is inconsistent with the distribution of 
rms values, we lower the velocity semiamplitude in the simulations 
and repeat the process until they are consistent. This particular 
example is for a trial period of 7.4 days with trial velocity semi- 
amplitude of 16.6 m s ~^ and 3 X 10* trials. The histogram shows 
the distribution of rms values derived from our Monte Carlo simu- 
lations while the dashed line marks the observed rms value of the 
residuals of HD 167042. The distribution of rms values migrates to 
smaller values as indicated by the arrow. This process continues 
for all relevant periods until the detection limits in period- velocity 
semiamplitude space are derived. 

C — B, and B — A) for each of the five target stars. The 
fowest values are close to 0.0010 mag, representing the 
standard deviations of the most stable comparison stars, 
which we take to be e. 

The photometric results for the five stars are summa- 
rized in Table [H Columns 4 and 5 give the duration of 
the photometric observations in days and the total num- 
ber of measurements, respectively. For all stars except 
6 Lyn, the duration of the photometric measurements ap- 
proximately equals or exceeds the orbital period of the 
planetary companion, though not by more than a couple 
of cycles. The standard deviations in columns 6 and 7, 
(T* and (Tc, vary somewhat from star to star due to differ- 
ences in factors such as stellar brightness, airmass, and 
seasonal photometric quality. Our method of computing 
the absolute variability of each target star takes these 
differences into account. 

Finally, the value of each star's absolute variability 
level, CT/, is given in column 8 of Table [H All values 
are significantly less than 0.001 mag (1 mmag) and are 
consistent with each star being constant. In cases where 
both the comparison stars and the target star are espe- 
cially stable, random errors in their measurements can 
result in < 0, i.e., cr/ becomes imaginary. This is the 
case for k CrB, for which we assume cr/ to be zero and 
the star to be constant. HD 210702 has the largest value 
of CT/ (0.00074 mag), though it is still less than 1 mmag. 
We classify HD 210702 as variable but append a colon to 
indicate some uncertainty. We find nothing significant in 
our periodogram analyses of all five stars. We conclude 
that there is no photometric evidence for brightness vari- 
ability levels in any of these five stars that could call into 
question the existence of their planetary companions. 



5. DETECTION LIMITS 

Information about the lack of planets around stars is 
critical for quantifying the occurrence rates of multiple 
Jovian planetary systems, computing frequencies of dif- 
ferent planet populations (jLagrange et al.l [2009a[) . and 
correcting for detectio n biases in radial velocity surveys 
(jCumming et al.ll2008[ ). We use Monte Carlo simulations 
of radial velocity measurements to assess the detection 
limits of stars in our sample. We derive upper limits for 
the velocity semiamplitude (K^p) as a function of period, 
which we convert to upper limits for minimum compan- 
ion mass as a function of semimajor axis using stellar 
masses. 

5.1. Method 

Our me thod for comp uting detection limits is similar 
to that of iLagrange et al . (2009a). We use the rms value 
of our radial velocities to determine which Keplerian or- 
bits of hypothetical planets are inconsistent with the ob- 
served rms values. For each star we begin with an ini- 
tial orbital period and a large velocity semiamplitude 
(-K'triai)- We generate a series of synthetic radial velocity 
measurements with random phases and circular orbits on 
the same dates that the observations were taken to cre- 
ate a distribution of rms values for that specific period 
and K trial- Noise is added to each synthetic measure- 
ment by drawing from a Gaussian distribution with a 
standard deviation equal to the total measurement-plus- 
jitter error at that date. If the distribution is inconsistent 
with the observed rms value at some confidence level (we 
chose 95.4% in this analysis), then iiTtriai is decreased 
and a new distribution is created (Figure [TT|) . If the ob- 
served rms value is consistent with the distribution then 
we adopt that iiTtriai value as the velocity semiamplitude 
upper limit, K^p, for that period. This process is then 
repeated over all relevant periods. 

5.2. Subgiants With Known Planets 

The detection of a planet around a star increases the 
probability that other planets are in that system, as find- 
ing a planet suggests that the stellar system in question 
is dynamically and was historically amenable to planet 
formation and rete ntion. Fo r know n planetary systems 
within 200 pc, Wri ght et all (|2009D find that additional 
planets are present at least 28% of the time. It is there- 
fore instructive to analyze the detection limits for the 
known planet-hosting stars in our sample so that certain 
regions of mass-semimajor axis phase space may be ruled 
out as not harboring additional planetary companions. 

To derive upper limits we generate 2500 artificial Ke- 
plerian orbits for each ii'triai- The if trial value is lowered 
in steps of 0.1 m s~^ until the rms distribution is con- 
sistent with the observed rms value at the 95.4% level. 
200 periods are evenly sampled from a log(P) distribu- 
tion between 2 to 5000 days. Results for the five IM 
subgiants are presented in Figure [12] and Table [5l 

Our observations are sensitive to planet minimum 
masses between a few times that of Neptune (0.053 Afjup) 
to several Jupiter masses for semimajor axes between 
0.04 AU to several AU. With a high probability we are 
able to exclude the existence of additional planets in 
the upper portions of each plot in Figure 1121 However, 
Neptune-mass planets in circular orbits cannot be ruled 



Retired A Stars and Their Companions. III. 



9 



TABLE 4 

Photometric Results for the Five Stars in Table [2] 
With Updated Orbital Parameters 







Date Range 


Duration 












Star 


APT 


(HJD - 2,400,000) 


(days) 


Nobs 


(mag) 


(mag) 


(mag) 


Variability 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


HD 167042 


12 


54128-55004 


876 


275 


0.00182 


0.00205 


0.00046 


Constant 


HD 192699 


12 


54192-55004 


812 


106 


0.00214 


0.00267 


0.00012 


Constant 


HD 210702 


12 


54100-54994 


894 


106 


0.00218 


0.00253 


0.00074 


Variable: 


ft CrB 


4 


49094-50250 


1156 


222 


0.00094 


0.00105 


0.00000 


Constant 


6 Lyn 


12 


54437-54928 


491 


238 


0.00153 


0.00163 


0.00011 


Constant 



10.0 i ■ ■ ; — i 10.0 




Fig. 12. — Detection limits for five planet-hosting subgiants based 
on Monte Carlo simulations of radial velocity measurements. The 
solid black and gray curves show the 95.4% and 68.3% minimum 
mass upper limit as a function of semimajor axis. Triangles show 
the locations of the known planets in these systems. The dashed 
lines indicate the semimajor axes that correspond to the baseline 
of the observations for that particular target. Planets located at 
larger semimajor axes will have completed less than one full period 
from the first to the last radial velocity measurement, hence the 
rise in the detection limit curves in that region. 

out by our observations. In general our detection lim- 
its are sensitive to the observed rms value, the internal 
measurement errors, and the assumed stellar jitter level. 
It is therefore unlikely that more observations will signif- 
icantly improve the resulting detection limits as stellar 
jitter is the dominant contribution to the rms budget. IM 
main sequence stars and giants both have higher mean 
jitter levels than subg iants (e.g.. iLagrange et al.ll2009at 
iSetiawan et"aI1 l2003bfi ; transit or microlensing surveys 
may therefore be the only tools available to study sub- 
Ncptune-mass planets in circular orbits around IM stars, 
although the larger radii and shorter evolutionary life- 
times of IM stars will make the detection of low-mass 
planets difficult even using these techniques. 

5.3. Subgiants from the Lick Survey 

Overall our sample contains a variety of radial ve- 
locity characteristics (see the "Notes" column of Ta- 
ble O for a summary): five targets are known planet- 
hosting stars, nine exhibit linear velocity trends, and 
four have large rms velocity scatter > 900 m s~^ (which 
roughly corresponds to the semiamplitude of a 13 Mjup 
object orbiting a 1.5 Mq star at 0.1 AU). The tar- 
gets with large velocity scatter are all known spectro- 
scopic binaries and include HD 8375 (discussed in §3.1), 



HD 65938 (Pourbai x et~an[20 04: Massaro tti et al.ll2008l). 
HD 179799 (.Setiawiaret al.ir2003b : .Setiawan et al.ll2004l : 
iMassarotti et a,l.l 120081) ."and HP 210211 (spectroscopic 
triple svstem: 'Horch et al.l l2002t iPourbaix erall 120041 : 

[Massarotti ct al. 2008). 

The two systems HD 90043 and HD 200964 both ex- 
hibit long-period (~ 1.5 yr) radial velocity variations 
which are likely caused by planetary-mass companions. 
Both stars are chromospherically quiet and do not show 
strong photometric variability.^^ LS periodograms con- 
firm strong periodicities for both systems with false 
alarm probabilities < 0.1%. Our radial velocity coverage 
is currently insufficient to obtain unique orbital solutions 
so we withhold a detailed discussion of these objects for a 
future publication. However, we note that single-planet 
orbital solutions fail to accurately reproduce the observed 
variations; instead, two-planet solutions produce signifi- 
cantly better fits. The inner planet of HD 90043 has a 
minimum mass of ~ 2.5 Mjup and is located at ^ 1.4 
AU (F - 1.2 yr), while the inner planet of HD 200964 
has a minimum mass of ~ 2.0 A/jup and is located at ~ 
1.7 AU {P ~ 1.7 yr). No additional stars in our sample 
exhibit periodicities with FAP < 0.1%. 

We derive detection limits for our full sample using 
the same parameters as in iJ5.2l for all stars with more 
than five radial velocity measurements. The results are 
presented in Table El We subtract the best-fitting linear 
trends for stars exhibiting constant accelerations. Ex- 
cluding the large rms velocity stars, the median (mini- 
mum) mass upper limits for our sample are ^ {0.21, 0.34, 
0.50, 0.64, 1.3} Mjup at semimajor axes of {0.1, 0.3, 0.6, 
1.0, 3.0} AU at the 95.4% confidence level. Our sample of 
subgiants is therefore typically sensitive to Saturn-mass 
planets out to 0.1 AU and Jupiter- mass planets out to 
^ 1 AU. Excluding stars with large rms velocity values, 
planets > 3 Afjup can be excluded out to 1 AU for all 
stars in our sample, and planets > 6 Mjup can be ex- 
cluded out to 3 AU for all but one of our stars. For the 
majority of our targets, however, our survey is sensitive 
to a few Jupiter-mass planets within 3 AU, or 1 Mjup 
planets within a few AU. 

6. JOVIAN PLANET FREQUENCY AROUND 
INTERMEDIATE-MASS STARS 

Several studies have provided hints that the occur- 
rence rate of Jovian p lanets increases with stellar mass. 
IJohnson et al.l (|2007af ) combine the results of the Califor- 
nia and Carnegie Planet Search around FGK and M-type 

Time-series photometric observations will be presented with 
the orbital solutions of HD 90043 and HD 200964 in a future pub- 
lication. 



10 



TABLE 5 
Detection Limits 







Other 






Mp sini (Afjup) for 


a 




rms 




Baseline 


Trend 




HD 


HIP 


Name 


(Mq) 


< 0.1 AU 


< 0.3 AU 


< 0.6 AU 


< 1 AU 


< 3 AU 


(m s"-"-) 


Nobs 


(yr) 


(m s-i yr-i) 


Notes^ 




Planet-Hosting Subgiants (Residuals) 


45410 


31039 


6 Lyn 


1.82 


0.15 


0.27 


0.37 


0.49 


0.86 


9.3 


64 


5.01 




P 


90043 


50887 


24 Sex 


1.91 


0.12 


0.20 


0.29 


0.37 


0.73 


8.4 


50 


4.37 




PP 


142091 


77655 


K CrB 


1.84 


0.07 


0.12 


0.17 


0.22 


0.37 


6.2 


62 


5.15 




P 


167042 


89047 




1.72 


0.11 


0.18 


0.25 


0.34 


0.59 


7.7 


68 


5.27 


2.1 


P, LT 


192699 


99894 




1.69 


0.20 


0.34 


0.49 


0.62 


1.40 


10.5 


50 


5.05 




P 


200964 


104202 




1.70 


0.13 


0.22 


0.33 


0.42 


0.81 


8.7 


91 


5.02 




PP 


210702 


109577 




1.85 


0.17 


0.32 


0.39 


0.51 


0.98 


8.8 


51 


4.76 




P 


Lick Observatory Subgiant Sample 


587 


840 




1.72 


0.20 


0.33 


0.42 


0.61 


1.05 


8.17 


22 


3.90 


-5.9 


LT 


8110 


6289 




1.64 


0.56 


1.04 


1.41 


2.13 


5.35 


21.3 


10 


2.02 






8375*^ 


6512 




1.66 


0.42 


0.68 


0.96 


1.25 


2.49 


19.0 


25 


4.05 


52.9 


SB, LT 


22682 


17049 




1.65 


0.34 


0.71 


0.88 


1.19 


1.98 


13.2 


10 


3.85 






25975 


19302 




1.59 


0.21 


0.32 


0.44 


0.63 


0.95 


7.71 


16 


4.69 






27536 


20263 


EK Eri 


1.87 


0.77 


1.36 


2.09 


2.95 


5.06 


32.6 


36 


3.24 






34538 


24679 




1.53 


0.28 


0.51 


0.70 


0.86 


1.66 


11.2 


12 


4.13 






37601 


26942 


24 Cam 


1.93 


0.28 


0.42 


0.68 


0.86 


1.54 


10.1 


17 


4.01 


9.2 


LT 


45506 


30815 




1.84 


0.16 


0.29 


0.43 


0.53 


0.84 


6.80 


14 


4.21 


-21 


LT 


57707 


35751 




1.81 


0.40 


0.91 


0.95 


1.17 


2.25 


10.8 


6 


2.24 






65938 


39198 




1.66 












2060 


4 


0.99 




SB 


73764 


42528 




1.91 


0.48 


0.73 


0.96 


1.33 


1.99 


13.1 


9 


3.30 






103484 


58110 


6 Vir 


1.98 


0.70 


1.10 


1.62 


2.71 


3.90 


25.8 


40 


4.09 






109272 


61296 




1.83 


0.22 


0.62 


0.70 


0.73 


1.26 


8.81 


16 


3.87 






111028 


62325 


33 Vir 


1.69 


0.19 


0.34 


0.47 


0.60 


1.07 


9.28 


23 


4.03 


-9.9 


LT 


123929 


69185 




1.61 


0.50 


0.59 


1.56 


2.40 


2.19 


12.3 


7 


2.62 






135944 


74690 




1.62 


0.17 


0.31 


0.43 


0.52 


0.95 


7.68 


16 


4.10 


53 


LT 


153226 


82989 




1.76 


0.19 


0.33 


0.50 


0.64 


1.04 


7.85 


18 


3.93 






179799 


94521 




1.74 


42.8 


85.0 


129 


177 


355 


2510 


8 


3.67 




SB 


184010 


96016 




1.80 


0.37 


0.64 


0.93 


1.13 


2.35 


15.2 


26 


4.07 






185351 


96459 




1.91 


0.17 


0.31 


0.43 


0.52 


1.11 


7.09 


17 


3.93 






202568 


104941 




1.57 


0.20 


0.31 


0.40 


0.54 


1.04 


7.10 


11 


3.82 


-5.2 


LT 


210211 


109281 




1.82 


125 


228 


330 


438 


794 


7930 


8 


2.25 




SB 


210404 


109338 




1.75 


0.68 


1.21 


1.77 


2.44 


14.8 


26.3 


11 


1.21 


-209 


LT 



Note. — Detection limits are for 95.4% confidence levels. Linear trends have been subtracted from targets with "LT.' 
"P" — planet; "LT" — linear trend; "pp" — probable planet, "SB" — spectroscopic binary (or higher order system). 



Detection limits are for the residuals after subtracting off the stellar companion and the linear velocity trend. 



stars with the frequency of planets using the entire Lick 
IM Subgiant Sample. They found that the frequency of 
Jovian planets appears to rise with stellar mass, reach- 
ing estimated rates of at least 9% for IM stars. Similarly, 
iLovis fc Mavorl ()2007n found a higher frequency of mas- 
sive planets (Mp sini > 5 Mjup) around IM stars than 
around solar-type stars. These previous estimates of Jo- 
vian planet frequency, however, have combined multiple 
published samples with different radial velocity sensitiv- 
ities and different sample selection criteria. While our 
sample size is considerably smaller in this work, we ben- 
efit from having a uniform data set observed over a long 
(~ 5-year) baseline with the same telescope and instru- 
ment setup. We also have the advantage of having well- 
characterized detection limits for the entire sample. 

Five stars in our sample are known planet-hosting sub- 
giants and we identify two more (HD 200964 and HD 
90043) as being strong candidates for having planetary 
companions (§5.3). The number of Jovian planet-hosting 
stars in our sample of 31 subgiants is therefore at least 
7. Bayes' Theorem can be used to derive the posterior 
pdf for the frequency of planets within ^ 3 AU: 



P(p|fc,n) (xP{k\p,n)P(j>), 



(5) 



where P(j)\k,n) is the posterior pdf for the probability 
of a star hosting a planet (p) given k detections in a 
sample of n stars, P(fc|p, n) is the pdf for observing k 
detections in the sample, and P{p) is the prior pdf for the 
probability that a star hosts a planet. The detection of 
a planet represents a Bernoulli trial, so P(k\p,n) follows 
a binomial distribution for the unknown parameter p. 
If we assume a uniform prior for P{p) between and 
1, the posterior pdf P{p\k,ri) simply follows a binomial 
distribution. 

The resulting posterior pdfs are shown in Figure 1131 
The median fraction of planet-hosting subgiants in our 
sample is 241^%, with upper and lower limits repre- 
senting the range encompassing 68.3% of the distribu- 
tion about the median. If we exclude subgiants exhibit- 
ing large rms velocity scatter, which is indicative of a 
stellar companion, then the fraction increases to 26lg%. 
These values can be compared to quantitative predictions 
from theoretical modeling of planet formation. For the 
mass rang e of our Lick su bgiant sample (1.5-2.0 Mq), 
[Kennedy fc KenvonI (|2008[ ) predict Jovian planet occur- 
rence rates to be ^ 10-15%. Our observed frequency is 
significantly higher, which may suggest that planet for- 
mation around IM stars is more efficient than previously 
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Jovian Planet Frequency 



N=28, k=7 
N=3^,k=7 



0.8 



1.0 



Fig. 13. — Probability density function (P(p)) for the frequency 
of Jovian-mass planets (p) around intermediate-mass stars interior 
to ~ 3 AU. Assuming a uniform prior, the posterior pdf for the 
probabiUty of harboring a planet follows a binomial distribution for 
k "successes" (planets detected) out of a sample size of targets. 
The dashed gray curve shows the distribution for the entire sample 
{N = 31 intermediate-mass subgiants), while the dashed gray line 
indicates the median of the distribution (24%). Several targets 
exhibited large radial velocity variations indicating a close stellar 
or brown dwarf companion; the black curve shows the posterior 
pdf excluding subgiants with large rms velocity scatter (28). The 
median value of this distribution is 26%. 

thought. The lower range of our results are, however, 
marginally consistent with the upper range predicted by 
[Kennedy fc KenyonI (|2008l ). 

Our results can be compared to the giant planet fre- 
quencies derived by other authors for solar-type stars. 
Our survey is sensitive to giant planets with m asses > 
1 .^'^jup within a few AU. John son et al.l (|2007a[ ) derive 
a planet frequency of ~ 4% for solar-type stars covering 
similar planet mass and semimajor axis ranges; the giant 
planet frequency they derive is significantly lower than 
the fre quency we find for IM stars. In a similar analysis, 
iLovis fc Mayor (2007) derive a frequency of ^^0.5% for 
planets with masses > 5 Mjup within 2.5 AU for solar- 
type stars. Over the same planet mass and semimajor 
axis range our program did not detect any planets. Using 
the binomial theorem, zero out of 28 stars translates into 
a 2 cr upper limit of 10.1% for the frequency of > 5 Mjup 
planets within 2.5 AU. Although poorly constrained, our 
frequency of the most massive giant pla nets around IM 
stars is consistent with the results from ILovis fc Mavod 
l|2007| ) for solar-type stars. 

Finally, we note that, intriguingly, the planet oc- 
currence rates for solar- and intermediate-mass stars 
resemble the occurrence rates of debris disks found 
throu gh excess infrared emission (^15% fo r Sun-like 
stars: iBryden et al.l[200a iTrilhng et a ll i20Q l >30% for 
IMstars: iRieke et al.ll2005l : iSu et alT l2006: M orales et all 
120091 ) ■ although at the moment no correlation has 
been foun d between the presence of planets and de- 
bris disks ("Moro- Martfn et al.ll2007l lK6spal et all 120091 : 
iBrvden et al.,.200a) . 

7. THE MASS-PERIOD DISTRIBUTION AND PLANET 
POPULATION SYNTHESIS 

The mass-period distribution of extrasolar planets 
around solar-type stars has been extensively studied 



in the literature (e.g., iTabachnik fc Trernaind 120021 : 
iLineweaver fc Gretheill2003l : iCumming et al.ll200a) . The 

simplest and most common parametric technique to 
model the distribution of giant planets is to fit a dou- 
ble power law to a mass and period histogram: 



dN oc A/"P'3rflnAfdlnP, 



(6) 



where dN is the number density of objects, M is the 
planet (minimum) mass, and P is the planet's orbital 
period. Values of a and /3 are gen erally consistent amon; 
authors; we take the values from lCumming et al.l (|200' 
as being representative of studies in the literature for 
Sun-like stars. Gumming et al. find a — -0.31 ± 0.2 and 
(3 ~ 0.26 ±0.1 for masses > 0.3 Mjup and periods < 
2000 days. Qualitatively this means that the mass and 
period distributions for planets around solar-type stars 
rise sharply (in linear space) for lower masses and shorter 
periods, respectively.^'* The probability of finding a low- 
mass, short-period planet is therefore much higher than 
finding a high-mass, long-period planet. The masses and 
periods of planets being discovered around IM stars is 
qualitatively quite different from those being discovered 
around solar-type stars. For example, out of > 20 planets 
discovered around IM stars, none orbit at semimajor axes 
less than 0.6 AU (~ 130 days for a 1.75 Mq star) and 
none have minimum masses below 0.70 Mjup (see Figure 
[Hand Table [TJ. This has generally been attributed to a 
different mass-period distribution for planets around IM 
stars, but the influence of higher jitter levels on detection 
limits in surveys targeting IM stars is usually not taken 
into account and no quantitative analysis has yet been 
performed. 

We use a Monte Carlo method to rigorously test the 
null hypothesis that the mass-period distribution for 
planets around solar-type and IM stars is the same. The 
question we seek to answer is the following: Given the 
values of a and (3 for solar-type stars, a planet occurrence 
rate, and the detection limits for our sample of 28 stars, 
what is the probability of finding 7 or more planets all 
having masses greater than 1.5 Mjup and semimajor axes 
greater than 1 AU?*^ 

We first test whether we can reproduce the observed 
number and properties of planets in our sample using 
the same Jovian planet frequency observed around solar- 
type stars. For each trial we randomly draw a new value 
of a and P from a Gaussian distribution centered on the 
values from Gumming et al. with a standard deviation 
equal to the quoted uncertainties. For each star in our 
sample we use the 10.5% planet frequency derived by 
Gumming et al. to determine whether that particular 

Equation|6]can also be expressed as dN oc P^^~^dMdP. 
In logarithmically-spaced bins the number of planets remains 
nearly constant. 

15 Three targets (HD 65938, HD 179799, and HD 210211) exhibit 
large rms velocity scatter and have either too few radial velocity 
measurements or too short time baselines to accurately fit orbit 
solutions, so the detection limits of the residuals are not available. 
We therefore exclude these targets from the simulation. 

In hypothesis testing, the probability of choosing extreme 
values (tail integrals) must be used rather than the probability of 
choosing a specific value. One could easily conceive of a discrete 
distribution in which the probability of choosing any particular 
value is small, so in this study we derive the probability of finding 7 
or more planets with masses greater than 1.5 A/jup and semimajor 
axes greater than 1 AU. 
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Fig. 14. — Comparison of the mass-period distributions of solar- 
type stars to intermediate-mass stars. Each histogram displays the 
results of a Monte Carlo planet population simulation Isased on 
values of a and f) from Cumming et al. ( 2008) following the power 
law distribution in Equation |6] an input planet frequency, and 
the detection limits from our sample. Using the planet frequency 
around solar-type stars (10.5%), the probability of detecting 7 or 
more planets is 0.024% (top left), while the probability of detecting 
> 7 planets with masses > 1.5 Afjup and semimajor axes > 1 AU 
is less than 10""^% (bottom left). When we adjust the input planet 
frequency so that each trial draws a new frequency value from a 
binomial distribution with 7 successes out of a sample of 31, 6.3% 
of the trials yield > 7 planets (top right), but only 0.002% of the 
trials yield > 7 planets with masses > 1.5 Afjup and semimajor 
axes > 1 AU (bottom right). 

star harbors a planet in our simulation. If a star harbors 
a planet, we randomly draw a mass and period for that 
planet from the power law distribution in Equation [51 
Finally, we check to see whether that planet would have 
been detected around that star based on the star's de- 
tection limits. We repeat this process 10^ times, saving 
the results after each trial. 

The results of our first Monte Carlo simulation are 
shown in the left panels of Figure [TH Out of 10^ tri- 
als, 24 produced 7 or more planets with any mass or 
semimajor axis, yielding a probability of 0.024% (Fig- 
ure [T3J top left). From the same simulation, however, 
no trials produced 7 or more planets with masses above 
1.5 Mjup and semimajor axes greater than 1 AU, yield- 
ing a probability of less than 10~^ (Figure [T4j bottom 
left). The input model values of a and [3 and/or the 
input planet frequency are therefore inconsistent with 
the number and properties of planets we observed in our 
sample of intermediate-mass stars. 

To determine whether increasing the planet frequency 
changes the results, we perform a similar Monte Carlo 
simulation except instead of the 10.5% occurrence rate 
we randomly draw a new planet frequency for each trial 
following the binomial distribution for 7 detections out 
of a sample of 28 (yielding a median planet frequency of 
~ 26%). For this exercise, the values of a and (3 are the 
same as the Cumming et al. values, but the probability 
of a star hosting a planet changes for each trial. The 
results of this simulation are presented in the right-hand 
panels of Figure [TD The probability of detecting 7 or 
more planets with any mass and semimajor axis rose to 
6.3%, but the probability of detecting 7 or more planets 
with masses above 1.5 Mjup and semimajor axes larger 



than 1 AU is 0.002%. Even when accounting for the 
different planet occurrence rates, we can rule out the 
Sun-like mass-period distribution for planets around IM 
stars at a confidence level of over 4 a. 

These simulations can easily be extended to exclude 
a wider range of a and [3 parameter space with varying 
degrees of confidence. We do this by running the same 
Monte Carlo simulation drawing from a binomial planet 
frequency distribution (with 7 detections out of 28 stars) 
over the following ranges of exponents: -2 < a < 8 (Aa 
= 0.5) and -2 < /3 < 8 (A/3 = 0.5). For each a-(3 pair 
we run 10'^ trials. The results are displayed in Figure 
[T51 each value of the grid represents the probability of 
detecting 7 or more planets with masses above 1.5 Mjup 
and semimajor axes greater than 1 AU. This technique 
allows us to exclude large regions of a- [3 space as being 
inconsistent with our observations. Negative values and 
small positive values of a and /3 are rejected with high 
confidence levels, while larger values of a and /3 are able 
to reproduce the observed properties more often. Quali- 
tatively this makes sense: we did not detect any low-mass 
(small a) or short-period (small (3) planets. We note that 
the inclusion of detection limits is critical in our analysis 
as it allows us to "detect" or "miss" simulated planets 
based on whether they fall above or below our detection 
threshold. 

8. DISCUSSION 

Little is known about the Saturn- or Neptune-mass 
population of planets around IM stars. Most planets or- 
biting IM stars discovered to date have masses > 2 Mjup 
(Table [T|), leading to suggestions that mo re massive stars 
tend to produce more massive planets (|Lovis fc MavoH 
I2007f ). It is unclear, however, whether this is a bona fide 
property or a result of an observational bias caused by 
higher jitter levels in IM stars which might mask the sig- 
nals of lower-mass planets. Moreover, few authors have 
assessed the limiting planet masses attainable by Doppler 
surveys of evolved IM stars. Despite the dearth of Jovian 
planets at small semimajor axes, could a population of 
Saturn-mass planets exist interior to 1 AU around IM 
stars? Would the higher jitter levels of evolved IM stars 
prevent the discovery of Neptune-mass planets? We ad- 
dress these questions in this study by deriving detection 
hmits for our sample of 31 IM (1.5 < M^/Mq < 2.0) 
subgiants. 

Typical detection limits for our stellar sample include 
planet masses down to - {0.21, 0.34, 0.50, 0.64, 1.3} 
Mjup within {0.1, 0.3, 0.6, 1.0, 3.0} AU at the 95.4% 
confidence level, excluding targets with large radial ve- 
locity scatter. We can therefore rule out the existence 
of hot Saturns within 0.1 AU and Jovian planets out 
to 1 AU for most stars in our sample. These detection 
limits suggest that the notably high masses of planets 
from our sample (> 2 Mjup) compared to planet masses 
around solar-type stars may be caused by a real differ- 
ence in planet population characteristics. For example, 
if a population of planets with random masses existed 
at 1 AU then we would expect an observational bias to 
result in an observed mass distribution that was trun- 
cated near the typical mass detection limit. As Jovian 
planets would have been detected out to ^ 1 AU, the 
higher planet masses uncovered so far may be indicative 
of a real trend. 
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Fig. 15. — Confidence regions for excluding pairs of a and /3 exponents for the parametric mass-period distribution dN oc M" d\nM d\nP 
based on planets in our sample. The values at each o-/3 pair indicate the fraction of the time our Monte Carlo simulations yield the high 
number (> 7) and the properties (Mp > 1.5 Mjup, a > 1 AU) of planets observed in our sample of intermediate-mass stars. We can 
exclude small values of a and /3 with a high level of confidence. The mass-period distribution for Sun-like stars from lCumming et al.l 1)2008 ) 
is labeled. 



We test this idea quantitatively by comparing the 
mass-period distribution of planets around Sun-like stars 
to the observed number and properties of planets in our 
sample. Even when correcting for the higher planet oc- 
currence rate found in our sample, the values of a and 
(3 in Equation [6] fail to reproduce the number (7/28), 
masses (> 1.5 Mjup), and semimajor axes (> 1 AU) 
of planets from our sample at a confidence level of > 
4 a. We conclude that the frequency and mass-period 
distribution of planets around IM stars is different from 
those around solar-type stars. Increasing the mass of the 
host star by a mere factor of 1.5-2 results in an entirely 
new planet population which is characterized by a high 
frequency (~ 26%) of high-mass planets (Afpsini > 1.5 
Miup) at large semimajor axes (a > 1 AU). 

The detection limits of the residuals of the planet- 
hosting stars in our sample show the strength of sub- 
giant jitter levels and therefore trace the sensitivity lev- 
els attainable to Doppler surveys targeting subgiants. 
The residual rms velocities range from ~6-10 m s~^. 
The detection limits indicate that planetary compan- 
ions with minimum masses above {0.2, 0.4, 0.5, 0.7, 1.4} 
Mjup within {0.1, 0.3, 0.6, 1.0, 3.0} AU cannot exist in 
these systems. If these detection limits are representa- 
tive of the typical jitter levels of subgiants then the higher 
rms velocity values observed in other subgiants from our 
sample may suggest as-yet-unrecognized low-mass com- 
panions. Given our jitter-dominated detection limits, 
the prospects of discovering Neptune-mass planets {'^ 



0.053 Afjup) in circular orbits around IM stars using the 
Doppler technique is not encouraging. If the dominant 
source of jitter in subgiants is from p-mode oscillations, 
observing strategies that include longer integration times 
or repeated exposures over hour- lon g timescales may hel p 
to partially overcome this hurtle (see| O'toole et all20Ml . 
Our ongoing survey of intermediate-mass subgiants at 
Keck Observatory will address this issue by using the 
high RV precision achievable with the HIRES spectrom- 
eter and larger telescope aperture of Keck compared to 
the Lick 3m. 

We derive updated orbit solutions using new obser- 
vations for five of the known planet-hosting IM sub- 
gia nts in our sam p le, whi c h were originall y announced 
bv I Johnson et al.N 2007bD. iJohnson et all (|2008l ). and 
iSato et al.l ()2008bl ). Our results are in excellent agree- 
ment with those previously reported in the literature. 
Our parameter uncertainties are typically a factor of 
2 smaller than the published values as a result of longer 
baselines and more radial velocity measurements. We 
also note that there is no evidence for periodicity in the 
residuals of the previously-known planet-hosting stars. 
We acquired time-series photometric observations of the 
five known planet-hosting stars with the Automatic Pho- 
tometric Telescopes at Fairborn Observatory. We find no 
evidence for brightness variation levels in any of the five 
stars that could call the existence of their planetary com- 
panions into question. 

The eccentricities of three of the five planets in our 
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sample with accurate orbital solutions are consistent with 
zero and emphasize a low-eccentricity trend for plan- 
ets around IM subgiants, all but two of which have 
eccentricities below 0.3 (Table [ij. This is in contrast 
to planets around solar- type stars, which have approx- 
imately uniform eccentricity distributions between 0.0 
and 0.8 for semima jor axes > 0.3 AU (Butler ct al. 200^; 
I Wright et all 120091) . Interestinglv. IWright et al.. (■2009f ) 
find that the eccentricity distributions for planets with 
masses < 1 Mjup peaks at e < 0.2, while the planets 
with masses > 1 Mjup have more uniformly distributed 
eccentricities between 0.0 < e < 0.6. Planets around IM 
stars show the opposite trend: high-mass planets tend to 
have circular orbits. We will explore this effect in detail 
in a future publication in this series. 

Information about planets orbiting IM stars has been 
driven by observations rather than theory. There is a 
growing need for stellar mass-dependent theoretical mod- 
els of planet formation that make testable predictions 
about the physical and orbital characteristics of planets 
in this stellar mass regime. Specifically, the abundance of 
low-mass planets, the eccentricity distribution, and the 
fraction of multiple planetary systems are but a few of 
the many outstanding questions to be addressed in this 
young field. 

We also encourage high-contrast imaging campaigns 
to include more A-type stars in their surveys. Five plan- 
ets have already bee n directly imaged o rbiting A stars: 
HR 8799 b, c, and d (jMarois et a l. 20081, Formalhaut b 
(iKalas et alll2008D. and /3 Pic b (jLagrange et al...2009bl . 
Moreover . iMarois et al.l ()2008[ ) made their discovery after 



observing only a few early-type stars, in contrast to the 
many hundreds of late-type stars that have yielded null 
detections. These results, combined with the higher in- 
ner planet occurrence rate we measure, suggest that more 
planet-hunting imaging surveys in this mass regime will 
yield fruitful results. 
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